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Preface 



This review first describes the evidence that strongly suggests the existence of the 
metal-insulator transition (MIT) in a two-dimensional electron system in Si regardless 
of the amount of disorder. Extensive studies of the charge dynamics demonstrate that 
this transition is closely related to the glassy freezing of electrons as temperature 
r — >■ 0. Similarities to the behavior of three-dimensional materials raise the intriguing 
possibility that such correlated dynamics might be a universal feature of the MIT 
regardless of the dimensionality. 
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0.1 Introduction 



The discovery of many novel materials over the last couple of decades has revived 
interest in the metal-insulator transition (MIT), one of the longstanding, fundamental 
problems of condensed matter physics. Many such materials, including manganites, 
diluted magnetic semiconductors and high-temperature superconductors, are created 
by doping an insulating host and thus find themselves close to a conductor-insulator 
transition. There is substantial evidence that, in many cases, strong electronic correla- 
tions (Mott localization) and disorder (Anderson localization) both play an important 
role in this regime (see [Miranda and Dobrosavljevic 2005| for review). In general, the 
competition between the Coulomb repulsion, which favors a uniform distribution of 
electrons (Wigner crystallization), and disorder, which favors a random one, leads to 
the frustration in the system. This means that electrons are unable to satisfy these 
requirements simultaneously, resulting in the absence of a unique ground state and 
the emergence of a large number of metastable states or minima ("valleys") in the 
free energy landscape. Metastable states, which correspond to different charge con- 
figurations, are separated by high barriers, leading to slow dynamics, divergence of 
the equilibration time, and breaking of ergodicity. Such a system is called a Coulomb 
glass. However, there are many other types of glassy systems. A common denominator 
in all of them is their complex or "rugged" energy landscape (Fig. 10. The most 
extensively studied and best known glasses in condensed matter physics are proba- 
bly conventional spin glasses ( [Binder and Young, 1986[ ) , such as Cu:Mn. On the other 
hand. Coulomb glasses were first anticipated theoretically almost three decades ago 
pavies et al, 1982|[G"runewald et al, 1982|[Pollak and Ortufio, 1982||Davies et al, 1984 



[PoUak, 1984 1 in situations where electrons are strongly localized due to disorder, 
but experimental studies have been scarce. Some of the recent work on electron 
glasses in such strongly localized regime, away from the MIT, has been reviewed by 



Amir et al. (2011). In the opposite limit of well-delocalized electrons, one expects a 
single, well-defined ground state and the absence of glassiness. Not surprisingly, it is 
the behavior in the intermediate region, near the MIT, that has been most difficult to 
understand. 

In the presence of disorder, the local electron density undergoes strong spatial 
fiuctuations. Therefore, it is plausible to expect that, in the vicinity of the MIT, the 
local density in some areas may be higher than the average and correspond to the 
metallic state, whereas in the remaining regions, the local density may be lower than 
the average and the electrons are localized. In fact, it has been proposed that, in 
weakly doped Mott insulators near the MIT, the system will settle for a nanoscale 
phase separatior0 between a conductor and an insulator ( [Gor'kov and Sokol, 1987| 



Kivelson et al, 2d03{ [Dagotto, 2002| . New powerful experimental techniques, includ- 



ing imaging by scanning tunneling microscopy (Kohsaka et al., 20071, have indeed 
firmly established the existence of charge inhomogeneities in strongly correlated elec- 
tron systems, such as cuprates ( [Orenstein and Millis, 2000{ [Dagotto, 2005D . This ob- 
viously leads to the possibility for a myriad of competing charge configurations and 
the emergence of the associated glassy dynamics. Recent studies ( [Raicevic et al., 2007[ 

^Global phase separation is not possible as it would violate charge neutrality. 
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[Raicevic et ai, 2008| [Jelbert et al., 200"8{ |Raicevic et ai, 201 ![ ) of a lightly doped, in- 
sulating La2-2:Sr2:Cu04, x = 0.03, have found several clear signatures of glassy charge 
dynamics as temperature T — ;> 0, consistent with an underlying glassy ground state 
that results from Coulomb interactions. Further work is needed, however, to see how 
this glassy state evolves as the system approaches the transition from an insulator to 
a conductor. 

It is also desirable to extend studies of charge dynamics to other types of materials 
to explore a possible generality of glassy freezing in strongly correlated systems near the 
MIT ( [Dobrosavljevic et ai, 2003l|Pankov and Dobrosavljevic, 2005| [Miranda and Dobrosavljevic, 2005[ ). 
Many materials do exhibit complex behavior due to the existence of several com- 
peting ground states ( [Dagotto, 2005[ ). In fact, the frustration caused by the com- 
petition of interactions on different length scales may give rise to glassy dynamics 
even in the absence of disorder ( [Schmalian and Wolynes, 2000[ ), while even a small 
amount of disorder may favor glassiness over various static charge-ordered states 



(Pankov and Dobrosavljevic, 2005). Even though the emergence of glassiness thus ap- 
pears to be ubiquitous at low temperatures. Coulomb glasses and out-of-cquilibrium 
systems in general remain poorly understood. Experimentally, studies of charge dy- 
namics near the MIT in many materials are often complicated by the accompany- 
ing changes in magnetic or structural symmetry, as well as by the glassy freezing 
of spins. Doped semiconductors, such as Si:P and Si:B, are free from such com- 
plications. They have been used extensively to study the critical behavior near the 
MIT ( [Miranda and Dobrosavljevic, 2005| [Sarachik, 1995 ) but, in spite of some early 
hints of glassiness in the insulating regime ( Paalanen et al, 1983[ [Monroe et al, 1987[ 
Monroe, 1990 1 , charge dynamics was not studied further until recently ( Kar et ai, 2003[ 



Thorsm0lle and Armitage, 2010[ ), as discussed below. 



Two-dimensional electron systems (2DES) in semiconductors (Ando et al, 1982), 
such as Si, provide another relatively simple system for exploring the interplay of 
electronic correlations and disorder (for a review of disordered electronic systems, see 



Free energy landscape 



Glass 




Configuration space coordinate 



Fig. 0.1 The free-energy (F) landscape of a glassy system. The horizontal axis represents 
the one-dimensional projection of the configurational coordinates of the degrees of freedom. 
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ILee and Ramakrishnan 1985^ . They have an additional advantage that aU the relevant 
parameters, carrier concentration, disorder and interactions, can be varied relatively 
easily. This review will describe some of the work that has demonstrated that the 2DES 
in Si is an excellent model system not only for studying the MIT in two dimensions 
(2D), but also for investigating the nearly universal noncquilibrium behavior exhibited 
by a large class of both three-dimensional (3D) and 2D systems {e.g. spin glasses, 
supercooled liquids, granular films). In fact, the work on the 2DES in Si provides 
additional strong evidence that many such universal features are robust manifestations 
of glassiness, regardless of the dimensionality of the system. 

The dimensionality plays an important role in the MIT. In 2D, for example, 
the very existence of the metal and the MIT had been questioned for many years. 
Recently, considerable experimental evidence has become available in favor of such 
a transition. Some of that work has been described in several review papers {e.g. 
IKravchenko and Sarachik 2004[) with a focus on very clean (low-disordered) samples. 
This review will first extend that discussion to the cases of higher disorder in order 
to (a) demonstrate that the 2D MIT is observed in all 2DES in Si regardless of the 
amount of disorder, (b) point out important similarities to the MIT in 3D systems, and 
(c) set the stage for the discussion of glassy dynamics near the MIT. As summarized 
below, the 2DES in Si exhibits all the main manifestations of glassiness: slow, corre- 
lated dynamics; nonexponential relaxations; diverging equilibration time, as T — > 0; 
aging and memory. The results provide strong support to theoretical proposals de- 



scribing the 2D MIT as the melting of a Coulomb glass (Dobrosavljevic et al., 2003 



Thakur and Neilson, 1996 Thakur and Neilson, 1999{|Chakravarty et al., 1999 Pastor and Dobrosavljevic, 1999| 



Dalidovich and Dobrosavljevic, 2002 1. The review concludes by further comparison to 
other complex materials, both 2D and 3D, and by identifying some open directions for 
future research. 



0.2 Metal-insulator transition in two dimensions 

In 2D systems in semiconductor heterostructures, the electron density Us can be varied 
easily over two orders of magnitude simply by applying voltage Vg to the so-called gate 
electrode. At low n^, the 2DES is strongly correlated: rs 1, where Ts = Ec/Ep oc 

— 1/2 

Us {Ec is the average Coulomb energy per electron and Ep is the Fermi energy; see 
also IKravchenko and Sarachik 20041 for more details). At the same time, the electrons 
"feel" a random potential (disorder) caused by charged impurities that are located 
away from the 2DES. Some of the striking experimental evidence for the MIT that 
occurs in a variety of "clean" , i. e. low-disordered 2D systems has been presented also 
in the chapter by S. V. Kravchenko (this volume). 

This section reviews work that demonstrates the existence of the 2D MIT regardless 
of the amount or type of disorder. The disorder, however, does affect the values of Uc, 
the critical density for the MIT, and the precise form of the temperature dependence 
of conductivity cr(T). For the sake of clarity, it is important to recall that a qualitative 
distinction between a metal and an insulator exists only at T = 0: cr(T = 0) 7^ in 
the metal and a{T = 0) = in the insulator. The experiments were performed on 
a 2DES in Si metal-oxide-semiconductor field-effect transistors (MOSFETs), where 
the low density regime can be reached more easily compared to other semiconductors 
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Fig. 0.2 High-mobility sample {fi ~ 2.5 m^/Vs): resistivity p — 1/a vs. T for 
n^lO^^cm-^) =8.55, 8.70, 8.84, 8.99, 9.13, 9.27, 9.56, 9.71, 9.85, 9.99, 10.4, 11.2, 11.6, 12.9, 
14.0 (from the top) ( |Jaroszynski et al, 2002 &| . Insets: a schematic of the sample (top view), 
and activation energies vs. Us; Uc ~ nl (here ~ 17). 



( Kravchenko and Sarachik, 2004D . In Si MOSFETs, the (Drude) mobility /i = (j/{nse) 
peaks as a function of Ug because, at very high Ug that are not of interest here, the 
scattering due to the roughness of the Si-Si02 interface becomes increasingly important 
( |Ando et al., 1982] ). The peak mobility at 4.2 K is commonly used as a rough measure 
of the amount of disorder. 



0.2.1 Effects of disorder 

High-mobility (low-disordered) samples. In high- mobility Si MOSFETs where, roughly 
speaking, the 4.2 K peak fi > 2 m^/(Vs), the low-temperature resistivity p exhibits 
a strong, metallic drop {dp/dT > 0) with decreasing T for Us > n* (Fig. 10. 2p 
( [Kravchenko and Sarachik, 2004[ ). dp/dT changes sign at a density n* and becomes 
insulatinglike [dp/dT < 0) for Ug < n*. Even though dp/dT < does not necessar- 
ily imply that p diverges at zero temperature (i.e. that a{T = 0) = 0), the density 
n* has been often assumed to represent the critical density for the MIT. However, 
other, more appropriate methods have been also used to identify Uc. For example, ric 
was determined based on both a vanishing activation energy and a vanishing nonlin- 
earity of current-voltage characteristics when extrapolated from the insulating phase 



(Pudalov et al., 1993 Shashkin et al, 2001 ). It was established that rtc although 
a small but systematic difference of a few percent has been reported such that ric < n* 
( [Pudalov et al., 1998a[ [Altshuler et al., 2001[ ). Fig. 10.21 inset illustrates the use of ac- 
tivation energies, determined from the fits of the data at the lowest Ug and T to the 
simply activated fornH p cx exp{EA/kBT). 

In the vicinity of (or n*), the resistivity curves can be collapsed onto the same 
scaling function of a single parameter Tq, i.e. p (T, Ug) = pcf{T/T[)) ( [Kravchenko et al., 1994 ) 



^ Close enough to ric, the data can be fitted almost equally well to the variable-range hopping law. 
This gives the same values of Uc within the error of the fit. 
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To is the same function of (5„ = (n^ — nc)/nc on both the metaUic and the insulating 
sides of the transition, Tq oc \S„\^'^ {z - dynamical exponent, v — correlation length 
exponent), zv = 1.6 ( [Kravchenko et ai, 1995^ . Here Pc{nc) is the so-called "separa- 
trix" , the temperature independent critical resistivity curve. On the metallic side, the 
scaling does not work at higher T, in the regime where p{T) goes through a maximum 
and then decreases as T is raised further ( [Kravchenko and Sarachik, 2004[ ). On general 
grounds, the scaling of p(ns, T) at low enough T represents a strong indication for the 
existence of the metal-insulator transition at T = ( [Sachdev, 1999[ ). 

At somewhat higher in the metallic phase, it becomes apparent that, after the 



initial rapid drop, p{T) becomes a very weak function as T — >■ (Pudalov et al., 1998a 



Pudalov et ai, 1998 &p . In fact, it is so weak that, on a logarithmic scale, it appears to 



saturate. As a result of this "saturation" of p{T), single-parameter scaling fails at the 
lowest temperatures ([Pudalov et ai, 1998 a[ ). However, a careful analysis of the data 
shows ( [Kim et aZ.,lM8 ; Was hburn et ai, 1999a] [Washburn et al, 1999 6| that all of 
the a{ns,T) (or p{ns,T)) curves can be scaled according to the more general form 
( [Belitz and Kirkpatrick, 1994[ ) a{ns,T) = ac{T)f{T /Tq), where the critical conductiv- 



ity CTc = a{nc,T) oc T^, i.e. it vanishes as T ^ 0. In other words, the critical conduc- 
tivity is not the "separatrix" , but instead it belongs to the insulating family of curves: 
rir. < nl, consistent with other studies ( [Pudalov et ai, 1998a| [Altshuler et ai, 2001 ). 



The exponent x, however, is very small: x = 0.1—0.2 (Kim et al., 1998 Washburn et ai, 1999a 



Washburn et al, 19996 ), so that a direct observation of (Tc{T) is very difficult at ex- 
perimentally accessible temperatures. A small value of x is also the reason for the 
apparent success of the single-parameter scaling in the limited range of T and n^. 
It is important to note that scaling with x ^ 0, also observed in 3D materials 
near the MIT ( [Sarachik, 1995[ ), does not contradict ( [Belitz and Kirkpatrick, 1994[ ) 
any fundamental principle for 2D systems. Indeed, violations of "Wegner scaling" 



(Belitz and Kirkpatrick, 1994), where x = {D — 2)/z {D - dimensionality), were pre- 
dicted for certain microscopic models ( [Castellani et al., 1987| [Kirkpatrick and Belitz, 1994| 
[Belitz and Kirkpatrick, 1995[ ) with strong spin-dependent components of the Coulomb 
interactions. 

While the general scaling form describes the data satisfactorily as T — !• 0, it would 
be desirable to explore scaling and the precise form of crc{T) over a wider range of T. 
It turns out that the situation becomes simpler in samples with more disorder, where 
the onset of "saturation" of a (or p) in the metallic phase is pushed to higher T. 

Intermediate-mobility samples and effects of local magnetic moments. Although the 
metallic temperature dependence of conductivity, da/dT < {i.e. dp/dT > 0), is 
less pronounced in samples with moderate mobility (4.2 K peak p ^ 1 m^/Vs), it 
was demonstrated early on (Popovic et al, 1997) that scaling a{ns,T) = (Tcf{T/To) 



is obeyed with exactly the same value of the exponent zv = 1.6 as in high- mobility 
samples. While the value of Uc, defined as n*, was higher than in high-mobility deviceqj, 
the corresponding w 13 was still comparably large. 



study of Si MOSFETs with different peak mobilities indicated a systematic increase of ric with 
disorder; ric was defined as the "separatrix" n* l |Pudalov et al., 1998a[ |. 
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Fig. 0.3 Intermediate-mobility sample (/i ~ 1 m^/Vs; sample 12): conductivity a vs. T for 
different Us (top), and scaling of a with temperature for the same data (and other Us not 
shown) in the Us range given on the plot (bottom). S„ = (us —nc)/nc is the reduced density; 
ric « nl (the corresponding ~ 15). For > nj, the scaling works only for T > 1.1 K. The 
substrate (back-gate) bias on the sample was Vsub = -40 V dFeng et aZ.,T999| . 



By extending the measurements to lower T (Feng et al., 1999), however, the "sat- 
uration" of cr(T) was observed to take place in the metallic phase at T as high as 
~ 1 K (Fig. 10.31 top). Moreover, the separatrix acquired an insulatinglike {da/dT > 0) 
temperature dependence. As a result, the simple, single-parameter scaling around the 
separatrix works only at T > 1.1 K on the metallic side (Fig. 10.31 bottom). Thus 
the data arc qualitatively similar to those on high-mobility samples, indicating that 
Uc < n*. However, since (t(T) in the metallic phase (n^ > ric) docs not have a simple 
enough form over a sufficiently wide range of T, it is still difficult to make reliable 
extrapolations to T = 0. Fortunately, in Si MOSFETs it is also possible to change 
the type of the disorder in the same sample. As shown below, this results in a sur- 
prisingly simple, precise form of cr(T') over a very wide range of T, allowing reliable 
zero-temperature extrapolations and striking scaling behavior. 

For a 2DES in Si, the disorder is due to the oxide charge scattering (scattering 
by ionized impurities randomly distributed in Si02 within a few A of the interface) 
and to the roughness of the Si-Si02 interface ( |Ando et al., 1982] ). For a fixed n^, it is 
possible to change the disorder by applying bias Vsub to the Si substrate (back gate). In 
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particular, the reverse (negative) Vsub moves the electrons closer to the interface, which 
increases the disorder. It also increases the splitting between the subbands since the 
width of the triangular potential well at the interface is reduced by applying negative 
Vsub- Usually, only the lowest subband is occupied at low T, giving rise to the 2D 
behavior. In sufficiently disordered samples, however, the band tails associated with 
the upper subbands can be so long that some of their strongly localized states may be 
populated even at low Ug, and act as additional scattering centers for 2D electrons. In 
particular, since at least some of them must be singly populated due to a large on-site 
Coulomb repulsion (tens of meV), they may act as local magnetic moments. Clearly, 
the negative Vsub reduces this type of scattering by depopulating the upper subband. 

Therefore, by varying Vsub, it is possible to study the effect of local magnetic 
moments on the transport properties of the conduction electrons in a systematic and 
controlled way. It has been established ( |Feng et al., 1999[ ) that scattering by local 
magnetic moments suppresses the metallic behavior with da/dT < 0. Indeed, the 
data suggest that in the T limit, the da/dT < {i. e. dp/dT > 0) behavior is 
suppressed by an arbitrarily small amount of scattering of the conduction electrons 
by disorder- induced local moments ( |Feng et ai, 1999[ ) . Figures I0.3f top) and I0.4f a) 
illustrate the striking effect of local moments on cr(T) in the same sample: while the 
"usual" , metallic behavior with da/dT < is observed in the absence of local moments 
[Fig. 10.31 top]. a{T) curves become insulatinglike [da/dT > 0) for all Us after many 
local moments are introduced with Vsub [Fig. 10. 4f a)]. However, this does not necessarily 
indicate the destruction of the metallic phase. In disordered 3D metals, for example, 
it is well known that the derivative da/dT can be either negative or positive near the 
MIT ( |Lee and Ramakrishnan, 1985[ ). On the other hand, the metallic behavior where 
a decreases but does not go to zero (as expected for an insulator) when T — > is new 
and unexpected in 2D. 

The analysis of the insulatinglike a(T) curves in Fig.lEla) shows dFeng et al.,200T\ 
that they follow a very simple and precise form over a broad (two decades) range 
of T and n^: a{ns,T) = a{ns,T = 0) + A{ns)T'^ [Fig. E^Jb)]. The high quality 
of the fits allows a reliable extrapolation of a{ns,T ~ 0), whose finite (i. e. non- 
zero) values mean that, in spite of the decrease of a{ns,T) with decreasing T, the 
2D system is in the metallic state. In particular, the zero-temperature conductivity 
a{ns,T = 0) is a power law function of Sn [Fig. I0.4r c)]: a{ns,T = 0) oc (5^ (/i w 3), 
as expected in the vicinity of a quantum critical point ( |Goldenfeld, 1992[ ), such as the 
MIT. The power law holds over a very wide range of Sn (up to 5) similar to what 
has been observed (Rosenbaum et ai, 19801 in Si:P near the MIT. In addition, even 
though the MIT occurs at different Uc in different samples, the critical exponents fi 
are the same [Fig. 10.4( c)]. as expected from general arguments ( Goldenfeld, 1992[ ). 
It has been also demonstrated ( [Feng et ai, 200 fj ) that, near the MIT, the data obey 
dynamical scaling a{ns,T) = ac{T)f{T/5^'^) [Fig. I0.4r d)] with a temperature de- 
pendent critical conductivity CTc = a{ns = nc,T) oc T^ {zv = 1.3 ± 0.1, x ~ 2.6, 
/i = x{zi') = 3.4 ± 0.4), both in agreement with theoretical expectations near a quan- 
tum phase transition ( [Belitz and Kirkpatrick, 1994 1 and consistent with the extrapo- 
lations of cr(r) to T = 0. 
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Fig. 0.4 (a)-(d) Adapted from Feng et al. (200 1| ; the same intermediate-mobility sample 



(sample 12) as in Fig. |a3 but here Vsub = +1 V. (a) a{T) for 0.3 < n^lO^cm"^) < 3.0 
(bottom to top) in steps of 0.1 x lO^^cm^^. (b) The same a{T) data plotted vs. the 
lowest Ua = 0.7 X lO'^^cm"^ and 0.3 < T < 2.2 K. Other measurements show that this o"(T) 
holds at least down to 0.020 K ( |Feng et al, 2001|[Eng"ei al, 2002| ). (c) a{ns,T = 0) vs. 5„ for 
samples 12 (dots) and 9 (squares). The dashed lines are fits with the slopes equal to the critical 
exponent fi. At the MIT, the corresponding ~ 22 and 17 for the two samples, respectively, 
(d) Scaling parameter To as a function of \S„\ for sample 12; open symbols: Us < n^, closed 
symbols: Us > ric. The dashed lines are fits with slopes 1.4 ±0.1 and 1.32 ±0.01, respectively. 
Inset: scaling of raw data ct/ctc ~ a/T^ in units of e^//iK^'^® for all Ua shown in (a) and 



T < 2 K. (e) Adapted from |Sjostrand and Stiles (1975[ l and |Sjostrand et al. (1976[ ): (t{T) vs. 

for three different samples and densities. The 4.2 K peak mobility was between 0.15 and 
0.33 mVVs. 
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The form of cr(T) is well established for metals containing local magnetic mo- 
ments, and it is believed to result from the Kondo effect ( [Hewson, 1993[ ). In fact, there 
is no other known mechanism that results in an increase of a as T^. Here this fea- 
ture provides the most direct evidence of the presence of local magnetic moments. It 
should be noted that, in general, one expects the behavior for a quantum impurity 
embedded in a Fermi liquid in any dimension. While the nature of this novel metallic 
state in 2D may require further study, its simple (t(T) allows for an unambiguous 
extrapolation to T = 0. The zero-temperature conductivity cr(ns,T ~ 0) decreases 
continuously, and follows a distinct power-law behavior as the MIT is approached. 
In particular, metallic a as small as lQ~^e? /h has been observed ( [Feng et ai, 200 
in a striking contrast to anything that has been reported in other 2D systems when 
da/dT < 0. A similar observation in 3D systems ( [Rosenbaum et ai, 1980[ ) has demon- 
strated the absence of minimum metallic conductivity, and has had a profound impact 
on shaping the theoretical ideas about the MIT. 

2DES in Si MOSFETs have been studied extensively for more more than four 
decades, so it is interesting that the behavior has been identified only relatively re- 
cently ( [Feng et ai, 200 f| ). A thorough examination of the early literature reveals, how- 
ever, that the samples discussed here are representative of a broad class of Si MOSFETs 



historically (and somewhat unfairly) known as "nonideal" samples ( Ando et ai, 1982). 
"Nonideal" samples could be made more "ideal" by applying Vsub and vice versa 
dAndo et al., 1982"] ), consistent with recent studies ( |Feng et al., 1999[|Feng et ai, 200l[ ) 
discussed above. Figure 10.4( e) shows an example of a{T) measured on samples with 
a modest peak mobility between 0.15 and 0.33 m^/Vs ( [Sjostrand and Stiles, 1975] 
[Sjostrand et ai, 1976[ ). The "saturation" (on a log scale) of cr{T) observed for T < 1 K 
was puzzling at the time, but the replotted data show clear behavior. 

Low-mobility (highly disordered) samples. Samples with very low 4.2 K peak mo- 
bility (fi < 0.1 m^/Vs) have attracted less attention because they do not exhibit a 
pronounced, if any, da/dT < metallic behavior. However, they not only exhibit the 
2D MIT, but also other similarities to the behavior of"clean" 2DES, such as the onset 
of glassy charge dynamics (Scction l0.3p . In some ways, it is even advantageous to inves- 
tigate samples with a lot of disorder. For example, the relevant electron densities, such 
as ric, are pushed to higher values (higher Ep), so that it is possible to reach much 
lower effective temperatures T/Tp {Tp - Fermi temperature) than in high-mobility 
samples. Therefore, comparative studies of samples with varying amounts of disorder 
should provide valuable insights into the physics of systems near the MIT. 

Detailed studies of transport and electron dynamics near the MIT have been per- 
formed on a set of Si MOSFETs with the 4.2 K peak mobility of only 0.06 m^Vs 
with the applied Vsub = — 2 V ( |Bogdanovich and Popovic, 2002 6[ ). This value of Vsub 
maximizes the peak mobility by removing the contribution of scattering by local mag- 
netic moments (see Sec. 10.2.11 above), at least in the experimental T-range. Because 
of the glassy fluctuations of a with time t at low and T (see Sec. 10. 3|) . the carrier 
density had to be changed at relatively high T (here ^ 0.8 K) in small steps in or- 
der to obtain reproducible values of the time-averaged conductivity (cr) . The behavior 
of {<j{ns,T)) (Fig. 10.51 left) turns out to be quite similar to that of high-mobility Si 
MOSFETs. At the highest rig, for example, the devices exhibit metallic behavior with 
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Fig. 0.5 Low-mobility sample {fi ~ 0.06 m'^/Vs) ( |Bogdanovich and Popovic, 2002^ . Left: 
time-averaged (a) vs. T for different Us- The error bars show the size of the fluctuations with 
time, rig, Ug, and ric are marked by arrows (rig - glass transition density). Right: {a) vs. T^'^ 
for a few ns (10^^cm~^) near iic. The solid lines are fits; the dashed line is a guide to the eye, 
clearly showing insulating behavior [{(t{T — > 0)) = 0]. 

d{a)/dT < 0. The change of (a) in a given T range, however, is small (only 6% for the 
highest Us = 20.2 x lO^^cm"^) as observed in other Si MOSFETs with a large amount 



of disorder ([Pudalov et ai, 1998a). Even though the density n* = 12.9 x 10^^ cm ^ at 



the separatrix (Bogdanovich and Popovic, 20026 Bogdanovich and Popovic, 2002a) is 
much higher than in "clean" samples, the value of (cr(n*)) = 0.5 /h is similar, which, 
according to Drude formula, corresponds to kpl < 1 {kp - Fermi wave vector, I ~ mean 
free path). Likewise, at the lowest Hs, the data are best described by the simply acti- 
vated form {a) oc exp^—EA/kBT). 

Perhaps the most striking difference between high- and low-mobility samples first 
becomes apparent when the vanishing of activation energy Ea is used to determine 
Tic- While in "clean" samples this gives Uc ^ n* (Sec. I0.2.ip . here Ea vanishes at 



ric ~ 5 X 10 cm~ , which is more than a factor of two smaller than n* (Fig. 10.5 
left). This suggests that there is a wide range of Ug on the metallic side of the MIT 
where o'(T) is insulatinglike. Indeed, close to ric, the best phenomenological fit to the 
data is the metallic power-law behavior (cr(ns,T)) = a{ns) + b{ns)T^ with x ~ 1.5 
(Fig. 10.51 right) ( [Bogdanovich and Popovic, 20026] ). The fitting parameter a{ns) is rel- 
atively small and, in fact, vanishes for ns(10^^cm~^) = 4.72 and 4.92. Therefore, 
Uc = (5.0 ± 0.3) X lO^^cm^^ (r^ ^ 7) based on the data on both metallic and insulat- 
ing sides of the MIT. Of course, a simple power law {a{nc,T)) oc is consistent not 
only with general expectations near the MIT and the behavior observed in 3D sys- 
tems ( [Belitz and Kirkpatrick, 1994 ), but also with the careful analysis of high- mobility 
2DES dKim et ai, 1998| [Washburn al, 1 999 a| [Washburn et ai, 19996D (Sec lOXT]) 
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and those with local magnetic moments ( [Feng et ai, 2001] ) fSec. l0.2."T]) . Here the expo- 
nent X takes a distinctly different value, presumably reflecting the different universality 
classes of those situations. 

The surprising non- Fermi liquid T"^/^ behavior is consistent with theory ( [Dahdovich and Dobrosavljevic, 2002 ) 



for the transition region between a Fermi liquid and an (insulating) electron glass. 
Indeed (see Sec. 10. 3p . the transition into a charge (Coulomb) glass in low-mobility 
samples takes place as T — > at a density rig, such that ric < Ug < n* (Fig. 10.51 left). 
The r3/2 correction is characteristic of transport in the intermediate, < < n„ 
region where the dynamics is glassy, but where a is still metallic [a{T — !• 0) 7^ 0] albeit 
so small that kpl < 1. Such "bad metals" include a variety of strongly correlated 
materials with unusual properties (Emery and Kivelson, 19951. Interestingly, the T^/'^ 



behavior can be revealed also in high-mobility 2DES by applying a parallel magnetic 
field. 

0.2.2 2D metal-insulator transition in a parallel magnetic field 

Since magnetic field B applied parallel to the 2DES plane couples only to electrons' 
spins, it is often used to probe the importance of spin, as opposed to charge, de- 
grees of freedom. Some of the intriguing results that have been obtained in parallel 
B in the vicinity of the zero-field 2D MIT have been described in the review by 



Kravchenko and Sarachik (2004). One of the main issues has been the fate of the 
metallic phase in a parallel B. From the insulating side, the critical density ndB) 
can be determined by extrapolating to zero the activation energy and nonlinearity 
of current-voltage characteristics (Sec. I0.2TT|) . For Ug > nc{B), however, the metallic 
da/dT < behavior observed in high-mobility samples is suppressed by B, mak- 
ing it even more difficult to determine the critical density from the metallic side. 
Nevertheless, a careful analysis reveals exactly the same metallic T^/"^ correction in 
high-mobility samples in parallel magnetic fields as in highly disordered samples in 
zero magnetic field, confirming the existence of the MIT in those two cases. 

Figure RTCT a) shows the {ug, B,T = 0) phase diagram obtained for a high- mobility 
sample ( [Jaroszyhski et ai, 2004 6[ ). The critical densities nc{B) were first determined 
using the activation energy method. Good agreement was found with the nc{B) de- 
pendence obtained on very similar samples using both activation energies and nonlin- 



ear current- voltage characteristic^ (Shashkin et ai, 2001). At low fields, ric increases 



with B, and then it saturates for 5 > 4 T, consistent with other studies that show 
that the 2DES is here fully spin polarized ( [Okamoto et al., 1999[ [Vitkalov et ai, 2000[ 
[Tutuc et ai, 2001D . Even though a{T) is very weak at higher n^, it is interesting 
to attempt to determine the separatrix, where da/dT = 0. The corresponding den- 
sity n*{B) > nc{B) and, most intriguingly, within the measurement error, it coin- 
cides with ng{B), the glass transition density ( [Jaroszyhski et ai, 2004a ). ng{B) was 



determined independently, based on the measurements of the fiuctuations in a{t) 
( [Jaroszyhski et al, 2004 6| (Sec. lUXT]) . 

The important question to address is the form of <j{T) in the narrow nc{B) < 
Us < ng{B) region. Here the data are best described by the metallic {{(t{T = 0)) ^ 0) 

^Because of the small sample to sample differences in the amount of disorder, the data from 
IShashkin et al. 2001l h ave been shifted up by 0.85 X lO^^cm to make the nc{B = 0) values coincide. 
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Fig. 0.6 High-mobility Si MOSFET. (a) T = phase diagram in a parallel magnetic 
field ( Jaroszyhski et al., 20046] ). The dashed lines guide the eye. The Uc values are from 
[Jaroszyhski et al. (20046 1 (dots) and |Shashkin et al. (200l] ) (triangles). The glass transition 



takes place at ng{B) > nc{B), giving rise to an intermediate, metallic glass phase. The den- 
sity at the separatrix nl, ~ rig within the error for all B. (b) a vs. T^'^ for several Us in 
the metallic glass phase for B = 2 T (n3(10^°cm^^) = 11.9,11.6,11.3,11.2,11.0,10.9,10.7 
from top; nc{B = 2 T) = 10.67 x 10^°cm~'^) for the same sample. Dashed lines are fits, 
(c) Zero-temperature conductivity a{T — 0) oc obtained from the fits shown in (b). (d) 
The positive magnetoresistance for the same sample ( [Jaroszyhski et al, 20046| ) and different 
densities ns(10^''cm~^), as shown. A strong increase in MR refiects a magnetic-field driven 
MIT at nc{B). 

power-law form (cr(n,.,B,T)) = {a{ns, B,T ^ 0))+b{ns, B)T^-^ [Fig. [Oljb)], similar to 
what was observed in the metallic glassy phase of highly disordered samples at B = 
(Sec. I0.2T1) . The extrapolated T = conductivities go to zero precisely at nc{B), in a 
power-law fashion {a{ns,B,T — 0)) cx 6!^ with fi « 1.5 [Fig. 10.6( c): <5„ = ns/nc{B)~l] 
( [Jaroszynski et al., 2004 &| that is in agreement with theoretical expectations near a 
quantum phase transition (Belitz and Kirkpatrick, 1994). Interestingly, there is some 
evidence (Fletcher et al., 2001) of similar behavior at B = with /.j ^ 1 — 1.5, obtained 



by extrapolating to T = the "saturation" of <7{T) in the da/dT < regime of different 
high-mobility Si MOSFETs. The striking power-law behavior shown in Fig. lO.Gf c) and 
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the remarkable agreement between nc{B) obtained from a{T) on both insulating and 
metallic sides of the MIT are strong evidence for the survival of the MIT and the 
metallic phase in parallel B. 

In low-disordered ("clean") samples at i? = 0, the critical density ric n* ~ Ug 
[Sec. l0.2.Tl Fig. l0.6r a)] and the intermediate, metallic glass phase is practically absent. 
A parallel B, however, increases its width, allowing the emergence of the character- 
istic T'^/^ correction to a. The increase of both Ug and ric, and the broadening of 
the metallic glass phase with B can be understood to result from the suppression of 
screening by a parallel B ( |Dolgopolov and Gold, 2000] [Herbut, 2001[ ), which mcreases 
the effective disorder. This, in turn, favors glassiness, consistent with theoretical ex- 
pectations (Dobrosavljevic et ai, 20031, and makes the behavior of "clean" samples 
more similar to that of highly disordered ones (Sec. I0.2T|) . The existence of the glass 
transition in high parallel -B, where the 2DES is spin polarized, provides evidence that 
charge, not spin, degrees of freedom are responsible for glassy ordering. This result 
clearly imposes a strong constraint on the types of theories that can be formulated to 
describe this phenomenon. Likewise, the broadening of the metallic glass phase with 
B also indicates that its existence is not due to spin. 

For rig near nc{B = 0), 2DES in various semiconductors exhibit a strong, pos- 
itive magnetoresistance (MR) ( [Kravchenko and Sarachik, 2004[ ), which has been a 
subject of great interest. The MR saturates at the field that corresponds to the 
full spin polarization ( |Okamoto et ai, 1999| [Vitkalov et al, 2000| [Tutuc et ai, 2001 



IVitkalov et ai, 2001D . As shown in Fig.ESJd) for the same high-mobility Si MOSFET 
discussed above, the MR jump is strong only for Ug not too far from nc{B = 0) since it 
takes place as the system undergoes a magnetic- field driven MIT at nc{B) [Fig. lO.Gf a)]. 
In the metallic phase where Ug ^ nc{B > 4 T), the MR is weak. 

In low-mobility Si MOSFETs, where the metallic glass phase is clearly observable 
already at i? = 0, the {ns,B,T = 0) phase diagram has not been studied. However, 
it is plausible to expect that the parallel B will broaden the metallic glass phase even 
further. 

Finally, in intermediate-mobility 2DES with local magnetic moments, it is relatively 
easy to map out the (ris, B,T = 0) phase diagram because of the simple form of cr(T) 
that holds over a wide range of T. In parallel B, a{ns,T) data (Fig. 10.71 top) are 
qualitatively similar to the B = case (Sec. I0.2.ip ( |Eng et ai, 2002[ ). At the lowest 
ns < Uc, for example, a decreases exponentially with decreasing T, indicating an 
insulating state at T = 0. For Ug > ric, o'(T) is weaker and its curvature is the opposite 
from the one expected for an insulating state. It clearly extrapolates to a finite value 
as T — )> 0, indicating a metallic phase (Sec. I0.2.ip . ric is identified as the density 
where CTc = c("-s ~ ricT) oc (Fig- 10.71 top), consistent with the B = case and 
in agreement with general arguments ( [Belitz and Kirkpatrick, 1994[ ). The exponent 
X = 2.7 ± 0.4 remains constant as a function of B. The critical density ric determined 
in this way at each given B allows the mapping of the {ns,B,T = 0) phase diagram 
(Fig. [QJl bottom). It should be noted that, at low fields (B < 2 T), a{ns,B,T) curves 



exhibit beautiful scaling (Eng et al, 2002) in agreement with T — > extrapolations 
and general considerations ( Belitz and Kirkpatrick, 1994] ) , providing additional strong 
evidence for a quantum phase transition in this system in parallel B. In analogy with 
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Fig. 0.7 Intermediate-mobility sample (/.i ~ 1 m'^/Vs) with local magnetic moments 
(V^sub = +1 V) Png et aL,l002l ). Top panel: cr(r) for sample 9 at B = 4 T. 
varies from 3.0 x lO^^cm"^ (top) to 0.7 x 10"cm"^ (bottom) in steps of 0.1 x 10"cm"^ 
nc{B = 4 T) = 1.4 X 10~^^cm~^ and is marked by the arrow. The dashed lines guide 
the eye. clearly follows a simple power-law dependence on T : (Jc oc T'^. Bottom panel: 
T = phase diagram for two samples. The dashed line guides the eye. The boundary 
between metallic and insulating phases is described by a power-law relation (see inset) 
[nc(-B) — nc(0)]/nc(0) oc (BjB*)^ at low fields, with the same crossover exponent /3 ~ 0.9 
for both samples (B* = 1 T for sample 9). Inset: the same data vs. B on a log-log scale. 
The dashed lines are fits with the slopes equal to /3. At B = 0, nc(10'^^cm-^) = 0.95 ± 0.05 
(ts ~ 17) and 0.85 ± 0.05 (r^ ~ 18) for samples 9 and 17, respectively. 

high-mobility samples, the MR exhibits a strong increase in the region of the magnetic- 
field driven MIT near ndB) ( |Eng et ai, 200T| ). For S > 2 T, where 2DES is spin 
polarized, nc{B) saturates, indicating that the MIT occurs between a spin-polarized 
metal and a spin-polarized insulator. The existence of the metallic phase at fields up to 
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18 T has been confirmed by the measurements of cr(T'), which retains a simple power- 
law form, albeit with a different exponent ( |Eng et al., 2002| . The charge dynamics has 
not been studied in these samples yet. 

It is interesting that the nc{B) dependence in samples with local moments (Fig. 
0.71 bottom) is quite similar to that obtained on high-mobility samples [Fig. I0.6f a)]. 
In both cases, at low fields ndB) increases with B in a power-law fashion: [nc{B) — 
nc(0)]/nc(0) oc . The crossover exponent /3 « 0.9 — 1 in Fig. 10.71 (bottom panel) 
and /3 = 1.1 ± 0.1 in Fig. I0.6f a). In other words, rig increases approximately linearly 
with B at low fields. It is remarkable that this dependence is essentially the same 
in both cases, even though the da/dT behaviors in the metallic phase at B = 
are strikingly different. The key features of the nc{B) phase diagram have been re- 
produced theoretically based on a scenario of quantum melting of a Wigner crystal 
as the mechanism of the MIT in sufficiently clean samples ( [Camjayi et al., 2008[ ). In 
general, a power-law shift of Uc with B is expected to occur in the case of a true 



MIT (Belitz and Kirkpatrick, 19941, and has been observed in several 3D systems 



(Rosenbaum et al, 1989 Bogdanovich et al, 1997 Sarachik et al, 1998 Watanabe et al, 1999) 



0.3 Glassy freezing of electrons in two dimensions 

Understanding the dynamics of glasses and other systems out of equilibrium is one 
of the most challenging and rapidly evolving topics in condensed matter research (see 
IBarrat et al. 2003p . Since vastly different types of systems exhibit similar behavior, 
it is tempting to search for common organizing principles and unified theoretical ap- 
proaches. However, despite some progress, there are still no well established theoretical 
frameworks for treating nonequilibrium behavior, which "remains largely uncharted 
territory" ( [Committee on CMMP 2010, 2007[ ). For example, although glassy behavior 
may dominate the low-temperature properties of many complex materials near quan- 
tum phase transitions ( Miranda and Dobrosavljevic, 200"5| (Sec. l0.1|) . such as the MIT, 
quantum glasses are even less understood than their classical counterparts. Experimen- 
tal studies of charge or Coulomb glasses ( Davies et al, 1982} [Griinewald et al., 1982[ 
[PoUak and Ortuho, 1982| [Davies et al., 1984| [PoUak, 1984|), which a re of particular 
relevance to the MIT, have been relatively scarce ([Monroe et al., 1987 Ben-Chorin et al., 1993 



Ovadyahu and Pollak, 1997 Vaknin et al., 1998 Vaknin et al, 2000 Vaknin et al, 2002 



Orlyan chik and Ovadyahu, 2004[| Ovadyahu, 2006a ||Ovadyah u, 2006 &[ [Martinez- Arizala et al., 1997[ 
Martincz-Arizala et a/., 1998; Biclcj ec and Wu, 20 01; Hcrnand ~e^ al, 2003[[Grenet, 20031 



Lee et al, 2005i|Grenet et al, 2007t[Kar et al., 2OO3(|Thorsm0lle and Armitage, 2010 1 
and mostly limited to insulating systems far from the MIT. Recent observations 



([Bogdanovich and Popovic, 2002 6[ Bogdanovich and Popovic, 2002a Jaroszyhski et al, 2002 & 
Jaroszyn ski et al, 2002a|[Popovic et al, 2003 [IJarosz ynski et al, 20 04&tlJaroszy nski et al, 2004a| 
Jaroszyhski and Popovic, 2006 Jaroszyhski and Popovic, 20076[ Jaroszyhski and Popovic, 2007a[ ) 



of glassiness in a 2DES in Si MOSFETs near the MIT open up opportunities for ex- 
ploring glassy phenomena in this important regime over a wide range of all the relevant 
parameters. 

There are two basic ways to explore the dynamics of glassy systems. The first one 
is to measure the response of the system to some kind of a perturbation. In a spin 
glass, for example, this would typically involve a study of the relaxation of the mag- 
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netization following some combination of rapid cooling and a change in the applied 



magnetic field (see [Vincent (2007 1 for a pedagogical review) . The second one is to 
measure the fluctuations of an observable with time (i.e. noise), which provides infor- 
mation on correlationfU. In spin glasses, transport noise measurements were required 
in order to provide definitive information on the details of glassy ordering and dynam- 
ics ( [Weissman, 1993| . Both approaches have been used to probe the dynamics of the 
2DES in Si, focussing on the conductivity a as the variable most relevant to the MIT. 
Measurements were performed on both high- and low-mobility Si MOSFETs, which 
also differ substantially in their geometry, size, and many other fabrication details, 
spanning essentially the entire range of Si technology. Thus the emergence of glassy 
dynamics proves to be a universal phenomenon in Si inversion layers, at least in the 
absence of disorder-induced local magnetic moments. The effect of local moments on 
charge dynamics still remains to be investigated. 

0.3.1 Fluctuations of conductivity 

The experimental protocol for measuring the a{t) fluctuations is simple, in princi- 
ple: the measurement of cr at a given Vg or electron density is set up at high 
temperatures, the sample is then cooled to the measurement T, and a is measured 
as a function of time. At the end of the measurement, the sample is warmed up 
to a temperature that is so high that a subsequent cool-down to the same mea- 
surement T would result in a reproducible value of the time-averaged (o'(i)) within 
the experimental uncertainty. The carrier density is changed at a high tempera- 
ture, and the protocol is repeated for different and T. In practice, the measure- 



ment set-up itself is often somewhat complicated (Bogdanovich and Popovic, 2002 & 
[Jaroszyhski et ai, 20026 Popovic et ai, 2003 ) because it is important to rule out var- 
ious extraneous effects, such as the fluctuations of T, Vg, or contacts, as possible 
sources of the measured noise. However, there are now a number of standard methods 



to accomplish this (see, e.g., Scofield (1987), Verbruggen et al. (1989)) 



In both low- and high-mobility samples, a exhibits strong fluctuations with time 
at low Us and T. Figure [0?8l left shows the fluctuations of (ct — (o'))/(cr) in a low- 
mobility sa mpldl for several at T = 0.13 K. It is striking that, for the lowest 
Us, the fluctuation amplitude is of the order of 100%. In addition to rapid, high- 
frequency fluctuations, slow changes over periods of several hours are also evident. 
The probability density function (PDF) of the fluctuations illustrates how the sample 
explores the free energy landscape (FEL). For Ug < rig, the PDF is always not only 
non-Gaussian but also has a very complex structure that changes with time (Fig. 10.81 
right). The PDF broadens with the increasing sampling time t, as the system has more 
time to explore the FEL. However, it explores it so slowly that it remains nonergodic on 
experimental time scales. For Ug > Ug, on the other hand, the PDFs become perfectly 
Gaussian on much shorter time scales, suggesting that the system reaches equilibrium. 

^In equilibrium systems, the connection between spontaneous fluctuations of a variable and the 
response of such a varia ble to a small perturbation in its conjugated field is given by the fluctuation- 



dissipation relation. See Leuzzi and Nicuwenhuizcn (2008) for the review and discussion of thermo- 
dynamics of out-of-equilibrium systems. 

®(. . .) represents averaging over time intervals of, typically, several hours. 
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Fig. 0.8 Low-mobility sample ( [Bogdanovicli and Popovic 20025] and Fig. lO.Sp . Left: Relative 
fluctuations of a vs. time for different at T = 0.13 K. Different traces have been shifted for 
clarity; the lowest Us is at the bottom and the highest Us at the top. Right: The color map 
of the probability density function (PDF) of the conductance (G) fluctuations as a function 



of sampling time t for the Uc < ns (10 cm 



5.58 < n„ data shown on the left. 



The amplitude of the fluctuations decreases dramatically from ^100% to less than 1% 
with increasing either Ug (Fig. 10.81 left) or T. 

An even more dramatic density dependence of the noise at low T is revealed by the 
study of the power spectra >S'(/) (/-frequency) of the relative changes in the conduc- 
tivity (cr(i) — ((t))/((t). Most of the spectra were obtained in the / = (10"'' - IQ-^) Hz 
bandwidth, where they follow the well-known empirical law S oc 1//" ( [Hooge, 1976[ 
[Weissman, 1988[ ). At the highest Us, S{f) docs not depend on Ug (Fig. 10.91 left) but, 
as rig is reduced below rig, S increases enormously, by up to six orders of magni- 
tude at low /. Moreover, for a given Us < Ug, S{f) increases exponentially with 
decreasing T (Fig. 10.91 left, inset). The observed dS/dT < makes it possible to 
rule out various simple models of noise (see, e.g., Bogdanovich and Popovic (20026), 



Popovic et al. (2003) for discussion). The most striking feature of the data, however, is 
the sharp jump of the exponent a at Ug ~ Ug (Fig. l0.9l right). While a ~ 1 ("pure" 1// 
noise) for > rig, a « 1.8 below Ug, reflecting a sudden shift of the spectral weight 
towards lower frequencies. Similar large values of a have been observed in some spin 
glasses above the MIT ( Jaroszyiiski et al., 1998| [Neuttiens et al., 2000D , and in sub- 
micron wires ( [Wrobel et al., 1998 1 . Both the increase in the magnitude of the noise at 
low / and the jump in a reflect a sudden and dramatic slowing down of the electron 
dynamics at Ug, indicating glassy freezing. The onset of glassy dynamics on the metal- 
lic side of the MIT 
glass phase for Uc < n 



i.e. at Ug > Uc 



5 X 10 cm~ , implies the existence of the metallic 
< Ug, which is consistent with the predictions of the model of 
interacting electrons near a disorder-driven MIT ( [Dobrosavljevic et al., 2003] ). Since, 
in the glassy phase, a decreases with increasing T (Fig. 10.91 right, inset), the large 
values of a and the jump in a(ns) are observable only at relatively low T. 
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Fig. 0.9 Low-mobility sample ( Bogdanovich and Popovic (20026 1 and Fig. 10. 5|) . Left 



The normalized noise power S{f ~ 3.16 X lO-^Hz) vs. Hs at T = 0.13 K 
( [Bogdanovich and Popovic, 2002a| . Below Ug ~ 7.5 x lO^^cm"^, the noise increases expo- 
nentially with decreasing ris. Inset: S vs. T for three different tIs (lO^^cm"^) given on the 
plot. Right: At Us « rig, the exponent a exhibits a sharp jump from ~ 1 at high Us 
("pure" 1/f noise) to ~ 1.8 at low Hs ( [Bogdanovich and Popovic, 2002^ . Inset: a vs. T 
for two different ns(10^^cm~^) (3.64 - open symbols, 3.86 - solid symbols) in the glassy 
phase ( [Popovic et al., 2003[ |. 

Qualitatively the same behavior has been observed in the resistance (or conduc- 
tance) noise of high-mobility samples, except that the glassy freezing takes place at 
Hg ?a ric ( [Jaroszynski et ai, 2002"6| . Importantly, in both types of samples, the char- 
acter of the noise also changes with density: while at low Ug both the "shape" and 
the variance of the noise exhibit random, nonmonotonic changes with time, at high 
enough Ug the noise always "looks" the same. This is illustrated in Fig. lO.lOr a) for a 
high-mobility sample. The figure shows (p — {p))/5p, where 5p = {{p — {p))^y^^, in 
order to make it easier to compare the signals. A quantitative measure of the spectral 
wandering with time, such as that observed at low ris, is the so-called second spectrum 
82(12, f), which is a fourth-order noise statistic. S2 is the power spectrum of the fluc- 
tuations of 5'(/) with time ( [Weissman, 1993[ Weissman et al., 1992[ ), i. e. the Fourier 



transform of the autocorrelation function of the time series of S{f). In general, if 
the fluctuators {e.g. two- level systems) are not correlated, -5*2 (/2, /) is white (indepen- 
dent of /2) dWeissman et ai, 1992[[\\^issman, 1988[ [Weissman, 1993[ ) and equal to the 
square of the first spectrum. Such noise is called Gaussian. On the other hand, ^2 has 
a nonwhite character, 52 oc l//2~^, for interacting fluctuators (Weissman et ai, 1992 
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Fig. 0.10 (a) Resistance noise in a high-mobility (low-disordered) sample for sev- 
eral Us shown on the plot (adapted from |Jaroszyiiski et al. 2002 (p — {p))/5p is 
plotted {Sp^ =variance, p - resistivity) in order to make the change in the char- 
acter of the noise with ris more apparent. Different traces have been shifted verti- 
cally for clarity. Exponent 1 — /3, a measure of correlations, vs. Us for (b) high-mo- 
bility [ric ~ 9.7 X lO^^cm"'^ ([Jaroszyh ski et al., 2002&t ] and (c) low-mobility samples 
[ric ~ 5.0 X lO^^cm"^ ( |Bogdanovich and Popovic, 2002fe[ |]. 



[Weissman, 1988| [Weissman, 1993| ). Therefore, the deviations from Gaussianity, or the 
exponent (1 — /3), provide a direct probe of correlations between fluctuators. Indeed, 
5*2 has been an important tool in studies of other glasses. 

A detailed dependence of the exponent (1 — /?) on Ug has been determined for both 
high- and low- mobility samples (Figs. lOnUl b) and (c), respectively) ( [Jaroszyhski et al., 20026 1 



In both cases, 5*2 is white [(1 — /3) = 0] for Ug > rig, indicating that the observed 1// 
noise results from uncorrelated fluctuators. It is quite remarkable that 5*2 changes its 
character in a dramatic way exactly at Ug in both types of samples. For < ng, 5*2 
is strongly non-Gaussian, which demonstrates that the fluctuators are strongly cor- 
related and provides an unambiguous evidence for the onset of glassy dynamics in a 
2DES at Ug. 

In the studies of spin glasses, the scaling of ^2 with respect to / and /2 has been 
used dWeissman et al., 1992[ [Weissman, 1993[ ) to unravel the glassy dynamics and, in 
particular, to distinguish generalized models of interacting droplets or clusters from 
hierarchical pictures. In the former case, the low-/ noise comes from a smaller num- 
ber of large elements because they are slower, while the higher-/ noise comes from a 
larger number of smaller elements, which are faster. In this picture, which assumes 
compact droplets and short-range interactions between them, big elements are more 
likely to interact than small ones and, hence, non-Gaussian effects and 5*2 will be 
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stronger for lower /. S'2(/2,/), however, need to be compared for fixed /2//, «■ e. on 
time scales determined by the time scales of the fluctuations being measured, since 
spectra taken over a fixed time interval average the high-frequency data more than the 
low- frequency data. Therefore, in the interacting "droplet" model, 5*2 (/2, /) should be 
a decreasing function of / at constant f2/f- In the hierarchical picture, on the other 
hand, 5'2(/2,/) should be scale invariant: it should depend only on /2//, not on the 
scale / dWeissman et al., 1992| [Weissman, 1993D . Figures iQ-llf a) and (b) show that 
no systematic dependence of 5*2 on / is seen in either high- or low-mobility 2DES 
( [Jaroszyhski et ai, 2002 The observed scale invariance of 5*2 (/2, /) signals that the 
system wanders collectively between many metastable states related by a kinetic hier- 
archy. Metastable states correspond to the local minima or "valleys" in the FEL (Fig. 
O.ip . separated by barriers with a wide, hierarchical distribution of heights and, thus, 
relaxation times. Intervalley transitions, which are reconfigurations of a large number 
of electrons, thus lead to the observed strong, correlated, l//-type noise, remarkably 
similar to what was observed in spin glasses with a long-range correlation of spin con- 
figuration (Weissman et ai, 1992 [Weissman, 1993[ ). On the other hand, in systems 



where both long-range and short-range interactions are present, such as lightly doped 
La2-a;Srj;Cu04, S'2(/2,/) reveals clearly ( [Raicevic et ai, 2008[ [Raicevic et ai, 2011[ ) 
[Fig. 10. lif e)] the presence of some characteristic energy scale, indicating that hier- 
archical models are not applicable and suggesting instead the formation of a cluster 
glass state. This is exactly what is expected in Coulomb systems with competing 
interactions (Sec. 10. ip . 

In the resistance noise measurements, the glass transition in a 2DES is thus man- 
ifested by (a) a sudden and dramatic increase of S{f) and a jump of a from « 1 to 
w 1.8, indicating the slowing down of the dynamics, and (b) a change of the expo- 
nent (1 — /3) from a white (zero) to a nonwhite (nonzero) value, indicating an abrupt 
change to a correlated statistics, consistent with the hierarchical pictures of glassy 
dynamics ( [Binder and Young, 1986t . For high- mobility 2DES in Si MOSFETs, low-T 
noise measurements were performed also in parallel B ( [Jaroszyhski et ai, 2004"6| . By 
adopting the same criteria for the glass transition in B, it was possible to determine 
ng{B) shown in Fig. ICGf a) and to establish that charge, not spin, degrees of freedom 
are responsible for glassy ordering (see Sec. I0.2.2|) . Measurements in both 5 = and 
i? 7^ show that glassy behavior generally emerges before the electrons localize {i.e. 
Ug > ric), consistent with theory ( [Dobrosavljevic et ai, 2003] ). The glassy signatures 
in the noise become gradually stronger as T decreases {e.g. Fig. l0.9l insets'). suggesting 
that the glass transition takes place as T — )■ 0. Strong evidence for Tg = and further 
support for Ug as the glass transition density is provided by the measurements of the 
response of the 2DES to a strong perturbation. 

0.3.2 Glassy response 

In MOSFET structures, the easiest and the most obvious way to excite the system 
is a sudden change of Vg. This method has been applied to several electron glasses 
([Ben-Chorin et ai, 19931 [Ovadyahu and Pohak, 1997[[Vaknin et aL , 19981 [Vaknin et ai, 2000 



Vaknin et ai, 2002[[Orlyanchik and Ovadyahu, 2004| |Ovadyahu, 2006aVOvadya hu, 20066| 



Martinez-Arizala et ai, 1997[ [Martinez- Arizala et 0171998 ; Grcnct, 2003 ; Grcnct ai, 20071 ), 
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Fig. 0.11 Scaling of tlie second spectra 52 measured at frequency octaves / — {fL,2fL) 
for a given Us in (a) liigli-mobility (for ris < ?ic) and (b) low-mobility (for ric < Us < ng) 
2DES ( |Jaroszyiiski et al, 2002fct . There is no characteristic time (or energy) scale observed, in 
agreement with hierarchical models. This is in contrast to (c) the insulating Lai.97Sro.o3Cu04, 
where S2 clearly shows the presence of some characteristic energy scale in that system 
( |Raicevic et al, 2008| |Raicevic et al., 2011[ ), suggesting that charge carriers form a cluster 
glass state. Dashed lines are linear fits to guide the eye. 



including a 2DES in Si ( Jaroszynski and Popovic, 2 006''Jarosz ynski and Popovic, 20076[ 
|Jaroszynski and Popovic, 2007a[|Jaroszynski and Popovic, 2009 ). So far, only low-mobility 
samples discussed above (Sees. [0?2T] and lO". 3.1 P have been probed in this way. In par- 
ticular, three different experiments have been carried out, as described below. 

Relaxations of conductivity after a rapid change of carrier density. A systematic 
study of the relaxations {i.e. time dependence) of the conductivity a{t) has been 
performed at different Us and T after a large, rapid change of rig (controlled by Vg) 
( [Jaroszynski and Popovic, 2006[ ). The sample was first cooled from high temperature 
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Fig. 0.12 Low-mobility sample with Ug ^ 7.5 x lO^^cm"^ and Uc « 4.5 x lO^^cm"^. 
(a) a{t) for V; = 11 V K(10"cm-2) = 20.26], V/ = 7.4 V K(10"cm-2) = 4.74], 
and T = 3.3 K. (b) Experimental protocol: Vg{t) and T{t). (c) All the data at 
short times, {i.e. before the minimum) are consistent with the Ogielski relaxation, 
a{t,T)/ao(T) oc (r;o„)"°(f/Tiot„)"° exp[-(t/r;o„)''] (dotted line). The lowest T data clearly 
deviate from the stretched exponential dependence (dashed line). The data have been col- 
lapsed with respect to the 2.4 K curve. Adapted from | Jaroszyhski and Popovic 2006[ 

(10 K) to the measurement temperature T with an initial gate voltage Vg. Then, 
at i = 0, the gate voltage was switched rapidly (within 1 s) to a final value Vj , 
and a{t,V/,T) was measured. In general, the results did not depend on any of the 
following: initial temperature, as long as it was > 10 K; Vg; the cooling time, which 
was varied between 30 minutes and 10 hours; the time the sample was kept at 10 K 
(from 5 minutes to 8 hours); the time the sample spent at the measurement T before 
Vg was changed (from 5 minutes to 8 hours) . Figures I0.12f a) and (b) show a typical 
experimental run with a(t > 0) exhibiting a rapid (< 10 s) initial drop followed by 
a slower relaxatior0. After going through a minimum, a(t) increases and approaches 
a value aQ{Vg^ ,T). A subsequent warm-up to 10 K and a cool down to the same 
measurement T = 3.3 K, while keeping the gate voltage fixed at Vg^ , shows that 
ao{Vj ,T) represents the equilibrium conductivity corresponding to the given Vg^ and 

^The RC charging times of the device and the circuit were at most 10 ms. 
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T. It is interesting that, even though initially it drops to a value close to ao, a first goes 
away from equilibrium before it starts approaching ctq again. At the end of the run, 
the sample was warmed up to 10 K, gate voltage changed back to the same Vg, and 
the experiment was repeated at a different T for the same V/. The whole process was 
then repeated for different values of Vj , in order to map out the density dependence 
of various relaxation parameters. It was established ( [Jaroszyhski and Popovic, 2006[ ) 
that, for < Ug, the relaxations have the following properties. 

After a sufficiently long t, a relaxes exponentially to its equilibrium value ctq- The 
corresponding equilibration time obeys the simply activated form r^q oc cxp{Eact/T), 
Eact ~ 57 K. While the microscopic origin of the activation energy Eact is not known 
yet, the activation to an upper subband in Si MOSFETs or to Si-Si02 interface traps 
has been ruled out ( [Jaroszynski and Popovic, 2007fe[ ). However, regardless of the equi- 
libration mechanism, the important result is that r^q — > oo as T — > 0, so that, strictly 
speaking, the system cannot reach equilibrium only at T = 0. In other words, the glass 
transition takes place at Tg = 0. Remarkably, a Monte Carlo study of the 2D Coulomb 
glass model has also found ( [Grempel, 2004[ ) an exponential divergence of Teg, signaling 
a glass transition at Tg = 0. There are, however, some differences in the detail between 
the model and the experiment, indicating a need for further refinement of theory. 

For short enough t, a{t) obeys a noncxponential, Ogiclski form ( [Ogielski, 1985 1 
ait,T)/ao{T) oc t"" exp[-(Vrio^)'^] (0 < aK) < 0.4, 0.2 < ^(n«) < 0.45) [Fig. c)], 
which is a product of a power law and a stretched exponential function. Both types of 
relaxations are considered to be typical signatures of glassy behavior and reflect the 
existence of a broad distribution of relaxation times. The scaling parameter tio^ cx 
exp{Ea/kT) {Ea ~ 20 K) so that, as T — > 0, tiow oo and the relaxations at- 
tain a pure power-law form oc t~°' . The Ogielski form, the divergence of tjo^, and 
the resulting power law relaxation at Tg are consistent with the general scaling argu- 
ments ( jHohenberg and Halperin, 1977D near a continuous phase transition occurring 
at = 0. These results are very similar to scaling observed in spin glasscfH above Tg 
(Pappas et ai, 20031. In a 2DES, a — ?> as ng ng, providing further evidence for 
Ug as the glass transition density. 

Tiow exhibits a very pronounced and precise dependence on the density: r;o„, oc 
exp(7ny^) (7-a proportionality constant). Since l/r^ = Ep/U oc 77.^^ in 2D (Sec. l0.2p . 
this particular form of Tiow{n-s) provides strong evidence that the observed out-of- 
equilibrium behavior is dominated by the Coulomb interactions between 2D electrons. 

Perhaps the most peculiar finding is that the 2DES equilibrates only after it first 
goes farther away from equilibrium. While the precise mechanism for nonmonotonic re- 
laxation remains controversial, studies of other materials [spin glasses ( [Jonsson et al., 1999| 
|J6nsson and Takayama, 2005 1 , manganites ( [Levy et al, 200"2 ) , insulating granular met- 
als (Kurzweil and Frydman, 20071, biological systems (Nelson, 2003 1] and some theo- 
retical models ( Morita and Kaneko, 2005] |Miyashita et al., 2007 1 suggest that it may 
be a general feature of systems that are far from equilibrium. 

It should be noted that it is quite remarkable that this one experiment has provided 
so many important results. To put things into a perspective, for example, it has not 



*Tg is finite in spin glasses. 
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been possible to determine Tg in other electron glasses so far, and scaling in spin glasses 
( [Pappas et ai, 2003[ ) has been observed only relatively recently, in spite of many more 
years of study. 

Relaxations of conductivity after a waiting time protocol: aging and memory loss. A 
key characteristic of relaxing glassy systems is the loss of time translation invariance. 



reflected in aging effects (Struik, 1978 Hodge, 1995 Rubi and Perez- Vicente, 1997). 



The system is said to exhibit aging if its response to an external excitation depends 
on the system history in addition to the time t. In a systematic study of the history 
dependence in a 2DES dJaroszynski and Popovic, 20076[ ), a{t) was measured after the 
temporary change of Us during the waiting time [Figs. lQTTSf a) and (b)]. The history 
was varied by changing ty^ and T for several initial (final) Ug. 

Two types of response have been observed: 1) monotonic, for relatively "small" 
excitations, where a{t) depend on t^, i.e. 2DES shows aging; 2) nonmonotonic, for 
sufficiently "large" excitations, where a{t) "overshoots" ctq (*-e- it first goes farther 
away from equilibrium) and relaxations no longer depend on tw (memory loss). The 
monotonic relaxations [Fig. lO.lST a)] are consistent with a power-law form at the short- 
est times (or lowest T) and, at the longest t, with a simple exponential approach to 
equilibrium. 

The aging is observed when <^ Tf.q{T), i.e. if the system is unable to equilibrate 
under the new conditions during [Fig. lO.l^ c)]. In that case, a{t) depends also on 
tyj: the system has a memory of the time t^. This is very similar to spin glasses, where 
T or B play a role analogous to that of n^. In the opposite case, when Teq{T) <C 
and the 2DES equilibrates at a new state, the relaxations do not depend on since 
the system, naturally, has no memory of the waiting time. Therefore, the overshooting 
is observed when the system is excited out of thermal equilibrium. This is analogous 
to the experiment described above (Sec. 10.3.2^ and thus sheds some light on that 
intriguing phenomenon. 

The gate voltage changes AVg employed in the relaxation experiments in a 2DES 
have been relatively large. For example, was changed up to a factor of 7, and thus 
the 2DES could go from the conducting to the insulating regime. Such large AVg 
are expected to trigger major rearrangements of the electron configuration since the 
corresponding shifts of the Fermi energjl^ AEp ^ ksT (Miiller and Lebanon, 2005). 



It is possible to speculate that such large perturbations might be somehow respon- 
sible for the peculiar overshooting effect. Considerable charge rearrangements, cou- 
pled with possibly substantial changes in the screening of the 2DES across the MIT 
( [Pastor and Dobrosavljevic, 1999[[Muller and loffe, 2004[|Pankov and Dobrosavljevic, 20"05l ), 
present a fundamentally different situation from the cooling of the 2DES at a fixed 
TLs when ksAT <^ Ep, where no relaxations have been observed. On the other 
hand, smaller perturbations of a Coulomb glass are expected to lead to memory 
effects ( [Lebanon and Miiller, 2005] ), in agreement with the observations for <C 
Teq{T). In that case, the final state has a large configurational similarity with the 
original state due to the shortness of t^. In InO^; films, another well-studied elec- 
tron glass, the overshooting of equilibrium has not been seen, but aging and mcm- 

^Ef[K]= 7.31r!,s[10"cm-2] for electrons in Si | |Ando et al, 1982^ . 
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Fig. 0.13 (a) a{t) for Vo = 7.2 V [ji^lO^cm^^) = 3.88 < Uc], Vi = 11 V 
[n^lO^^cm"^) = 20.26], t„ = 1000 s, and T = 3.0 K. (b) The corresponding experimen- 
tal protocol: Vg{t) and T(t). The results do not depend on the cooling time (varied from 
30 minutes to 10 hours), nor on ti and t2 (varied from 5 minutes to 8 hours each), (c) 
The symbols show the values of {T,tw) where the overshooting vanishes, i.e. the conditions 
that separate the two regimes. Open and solid symbols correspond to different samples; the 
blue dot-dashed line guides the eye. The purple dashed line represents Teq{T). Adapted from 



Jaroszyriski and Popovic (200761 



ory effects have been observed ( Ben-Chorin et al., 1993 Ovadyahu and Pollak, 1997| 
[Vaknin et ai, 1998{|Vaknin et ai, 2000 Vaknin et ai, 2002 ). Those experiments were 
done in the regime of small perturbations, because the typical change in the carrier 
density due to AV^ was ^ 1%, the system always remained deep in the insulating 
state, and T^q was satisfied. 

Aging effects across the MIT in 2D. Aging effects have been instrumental as a probe 
of complex nonequilibrium dynamics in many types of materials. In a 2DES, where 
the onset of glassy dynamics takes place on the metallic side, a study of aging, espe- 
cially across the MIT, is of great interest. Aging is observed if tw ^ Teq [Fig. 10.13( c)]. 
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Since the equilibration time T^q diverges exponentially as T — (Sec. 10.3. 2p . strictly 
speaking, the system can reach equilibrium at all T > 0. However, even at T that 
are not too low {e.g. ^ 1 K), T^q exceeds easily not only the experimental time win- 
dow but also the age of the Universe ( [Jaroszyhski and Popovic, 2006 ). This makes 
it relatively easy to study the out-of-equilibrium relaxation of a at times t <C Teq, 
where one expects to find properties common to other types of glasses. In strongly 
localized systems, such as InOx films, the aging function a{t,tw) is just a func- 
tion oft/tu, dBen-Chorin et al, 19931 [Ovadyahu and Pollak, 1997|[Vaknin et al, l"998l 
|Vaknin et ai, 2d00{ Vaknin et al, 2002 1. This is known as simple, or full aging. It is 
interesting that, in spin glasses, full aging has been demonstrated only relatively re- 
cently ( Rodriguez et al., 2003[ ). In general, however, the existence of a characteristic 
time scale tw does not necessarily imply simple t/tu, scaling ( [Barrat et al, 2003] ). In 
the mean-field models of glasses, for example, two different cases are distinguished: 
one, where full aging is expected, and the other, where no t/t^u scaling is expected 



(Bouchaud et al., 19971. Experimentally, departures from full aging arc common (Struik, 1978 



Rubi and Perez- Vicente, 19971. 



Aging was investigated in detail in a 2DES ( [Jaroszyhski and Popovic, 2007a| ) using 
the waiting time protocol [Sec. I0.3.'2l Figs. l0.14F a). (b)], but T was kept fixed at 1 K 
such that Teq was astronomical and, hence, the 2DES was always deep in the t^, r^q 
limit [Fig. I0.13f c')]. a{t,t.uj) were then explored systematically both as a function of 
final Ug and of the difference in densities during and after tw Figure [0?T4l c) illustrates 
the significant effect of t^ on (j{t). In fact, all the a{t,tw) data can be coUapsed onto 
a single curve simply by rescaling the time axis by tw [Fig. I0.14r d)]. Therefore, in this 
case, the system exhibits full aging at least up to < « (2-3)tto. The relaxations can be 
described by a power law a(t)/ao oc {t/tw)~°' for times up to about t^,, followed by a 
slower relaxation at longer t. This means that the memory of t^ is imprinted on the 
form of each a{t). 

For all Us < ric, a{t,tm) exhibit simple or full aging. However, as increases 
above ric, there is an increasingly strong departure from full aging. In some other 
glasses dStruik, 1978] [Alba et aZ.,l986l |Alba et aL,T987| , it was found that the data 
could be scaled with a modified waiting time (tu,)^, where is a fitting parameter 
(/i = 1 for full aging). Even though fi may not have a clear physical meaning, the 
//-scaling approach has proved to be a useful tool for studying departures from full 
aging ( [Struik, 1978| [Barrat et al., 2003t . By adopting a similar method in the study 
of aging in a 2DES, it was possible to achieve an approximate collapse of the data 
(Jaroszyhski and Popovic, 2007a). The plot of /i vs. uq [Fig. I0.14f e')] shows a clear 
distinction between the full aging regime for Ug < Uc, and the aging regime where sig- 
nificant departures from full scaling are seen. It is striking that the largest departure oc- 
curs at rig K, Ug. For rig > ng, some small relaxations are observed (only for kpl < 1) re- 
sulting from finite-temperature effects, but they vanish in the T — > limit (Sec. I0.3!2|) . 
It was also determined that /i does not depend on T ( [Jaroszyhski and Popovic, 2009[ ). 

These results are striking because they reveal an abrupt change in the nature of the 
glassy phase exactly at the 2D MIT itself, before glassiness disappears completely at 
a higher density Ug. In other words, this is strong evidence that the insulating glassy 
phase and the metallic glassy phase are different. Therefore, the difference in the aging 
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Fig. 0.14 (a) Vg{t) and r(i) in a typical experimental protocol, which always starts with the 
2DES in equilibrium at 10 K ( [Jaroszyhski and Popovic, 20076| ). (b) The corresponding o"(t). 
The relaxation of a during is too small to be seen on this scale, (c) a{t > 0) for several tu,, 
as shown; Vo = 7.0 V [no(10"cm-2) = 3.02 < Uc], Vi = 11 V [ni(10^^cm-2) = 20.26]. The 
dotted lines are linear fits for t < t^. (d) The same data as in (c) but plotted vs. t/tw The 
dotted line is a fit for t < with the slope —a — —0.081 ± 0.005. Inset: a{t = ls)/ao vs. t^. 
The dotted line is a fit with the slope a = 0.076 ± 0.005. (e) fj, vs. no for two samples, fi does 
not depend on ni. (f) Relaxation amplitudes vs. no for several t-w. Dotted lines guide the eye. 



Adapted from Jaroszyriski and Popovic (2007a I, and Jaroszyhski and Popovic (20091. 
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properties below and above Uc puts constraints on the theories of glassy freezing and 
its role in the physics of the 2D MIT. 

Furthermore, for a given tw, the amplitude of a{t)/ao has a peak at Us < ric 
[Fig. I0.14r f)]. reflecting an interesting and surprising suppression of the relaxations 
on the insulating side of the 2D MIT. While a clear understanding of this effect is 
lacking, it is plausible that collective charge rearrangements that are responsible for 
the slow dynamics will be suppressed as the 2DES becomes strongly localized. It 
would be also of interest to study aging deeper in the insulator, in the variable-range 
hopping regime, but that is not possible because of the small relaxation amplitudes and 
the large intrinsic sample noise. Of course, the effects of disorder could be explored 
further by extending the relaxation studies to cleaner (high- mobility) 2DES, where 
Ug > ric ( [Jaroszyiiski et ai, 2002 6| |Jaroszyhski et ai, 2004 . Regardless of the origin 
of the nonmonotonic behavior in Fig. 10.1 jT f). however, it is important to note that it 
reflects another change in the aging properties that occurs at the MIT. 



0.4 Summary 

Careful studies have provided strong evidence that all 2DES in Si MOSFETs exhibit 
a metal-insulator transition regardless of the amount or type of disorder (Fig. I0.15|) . 
This is confirmed by extrapolating (j{T) on both metallic and insulating sides of 
the MIT, and by dynamical scaling of a{ns,T) over a wide range of parameters in 
both B = and B 7^ 0. In B = 0, the scaling form cr(n^,r) = aciT)f{T/d^'') is 
obeyed with a temperature dependent critical conductivity Cc — (^{us — nc,T) (x T^, 
where the value of x depends on the disorder. Thus <Jc{T) belongs to the insulating 
family of curves. The results are consistent with general arguments near the MIT 
( [Belitz and Kirkpatrick, 1994 ) and they are similar to the behavior in doped semicon- 



ductors near the 3D MIT (Sarachik, 19951. The 2D metallic phase survives in very 



high parallel magnetic fields, long after the 2DES becomes fully spin polarized. 

Measurements of the charge dynamics in a 2DES have established that the glass 
transition takes place at Tg = for all < ng (Fig. 10.15]) . In general, the glassiness 
sets in on the metallic side of the MIT, i.e. at a density Ug > ric, thus giving rise to 
an intermediate, metallic glass phase between the "normal" metal (i.e. in the kpl > 1 
regime) and the glassy insulator. The glass transition and various glassy effects are 
observable only for kpl < 1, so that the intermediate phase is a very poor metal 
(cr(T = 0) 7^ <C e? /h). It is characterized by a very specific, non-Fermi liquid T"^/^ 
correction to a in both low-mobility samples in S = and high-mobility samples in a 
magnetic field. In high-mobility samples in i? = 0, the intermediate phase practically 
vanishes and the glass transition coincides with the MIT. 

The manifestations of the glass transition in a 2DES for < Ug, as demonstrated 
by resistance noise measurements, include a dramatic slowing down of the electron 
dynamics and correlated statistics consistent with the hierarchical picture of glasses. 
The results have been further confirmed and supported by the studies of relaxations, 
which provide evidence for the diverging equilibration time as T — >^ 0, broad distribu- 
tion of relaxation times in the system, and scaling (in the time domain), consistent 
with a continuous phase transition at Tg =0. The aging, one of the hallmarks of glassy 
dynamics, shows an abrupt change in its properties at the MIT, indicating that the 
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Fig. 0.15 Experimental phase diagram of the 2DES in Si MOSFETs. The metal-insulator 
transition takes place at the density ric and at T = in all samples regardless of the amount 
of disorder. The nature of the metallic phase at high Us, such that kpl > 1, is still under 
debate [see Kravchenko and Sarachik (2004 1]. The glass transition takes place at Tg = for 
all Ha < rig. In general, the glass transition is a precursor to the metal-insulator transition, 
I.e. Tic < Tig, giving rise to an intermediate metallic phase with a particular form of the 
non-Fermi liquid temperature dependence of conductivity. The aging properties of the glass 
change abruptly at ric, indicating different natures of the insulating and metallic glass phases. 
For sufficiently low disorder, the intermediate phase vanishes; Uc < ~ Ug . 

natures of the insulating glass and metallic glass phases are different. The 2DES in 
Si thus exhibits all the characteristics common to other out-of-cquilibrium systems, 
regardless of their dimensionality. Numerous similarities to spin glasses with long- 
range RKKY interaction arc particularly remarkable. Here, however, only long-range 
Coulomb interaction and potential scattering are present, and measurements in a par- 
allel magnetic field confirm that charge degrees of freedom, not spin, are responsible 
for glassy freezing. Therefore, the experiments show that the 2D MIT is closely related 
to the melting of this Coulomb glass. 



0.5 Discussion 

Experiments on a 2DES in Si clearly provide strong support to theoretical proposals 



describing the 2D MIT as the melting of a Coulomb glass ( Dobrosavljevic et ai, 2003| 



Thakur and Neilson, 1996 Thakur and Neilson, 1999{|Chakravarty et ai, 1999 Pastor and Dobrosavljevic, 1999 



Dalidovich and Dobrosavljevic, 2002 1. In particular, the model of the MIT as a Mott 



transition with disorder ( Dobrosavljevic et ai, 2003[ ) predicts the emergence of an in- 
termediate metallic glass phase in sufficiently disordered systems (Fig. l0.16"| . Changing 
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Fig. 0.16 Theoretical phase diagram (adapted from |Dobrosavljevic et al. 2003[ ). The Fermi 
energy and disorder on the two axes are expressed in units of the on-site interaction U . 
For large enough disorder, the metal-insulator transition is preceded by the glass transition, 
giving rise to an intermediate, metallic glass phase. Experimental trajectories, where the MIT 
is approached by varying Us, are represented schematically by the green curved arrows for 
the cases of low (trajectory 1) and high (trajectory 2) disorder. 

the carrier density in highly disordered samples, for example, corresponds tcF°l trajec- 
tory 2 in the theoretical phase diagram in Fig. 10.161 where the intermediate phase is 
predicted to be relatively broad. On the other hand, for low-disordered samples (trajec- 
tory 1), the metallic glass phase is very narrow or it vanishes, in agreement with exper- 
imental observations. The same theory predicts ( [Dalidovich and Dobrosavljevic, 2002[ ) 
the T^/^ correction to conductivity in the metallic glass phase, which is precisely what 
is found in the experiment. The emergence of the metallic glass phase, its dependence 
on the disorder, and the specific form of cr(T) in this regime were obtained using a 
mean-field theory approach, which is known to produce hierarchical dynamics and ag- 
ing in models of other glasses. Therefore, this approach seems promising in describing 
also the glassy charge dynamics in a 2DES in Si. There is currently no other theory 
available that predicts any of these experimental features. 

A molecular dynamics simulation of the crossover from a Wigner liquid to a Wigner 
glass in a 2D system of interacting electrons with disorder ( [Reichhardt and Reichhardt, 2004| ) 
has reproduced the main noise results (Sec. I0.3.1|) . In particular, a strong, orders of 
magnitude increase in the noise power and a jump in the exponent a was found at 
low T and rig, as well as the emergence of non-Gaussianity. By looking at the electron 
trajectories for a fixed period of time (Fig. I0.17| . it was found that, at the highest rig, 
electrons can flow freely throughout the sample. As ris is reduced, electron motion con- 
sists of a mixture of 2D and ID regions (Fig. lO.lTl upper right), but over longer times. 



^''Since Ep ~ Us and U ~ nj^^ in 2D and the disorder W ~ const, the change of Us in the 
experiment is described by (Ep/U) {W/U)~^ in the phase diagram. 
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Fig. 0.17 Electron trajectories for a fixed period of time for a fixed T and four electron 
densities, with tiie iiighest ni at upper left, 712 < ni upper right, 713 < n2 lower left, and the 
lowest n4 < 713 lower right. Adapted from Ref. ( [Reichhardt and Reichhardt, 200"4| ). 

the motion still occurs throughout the whole sample. At even lower Us (Fig. lO.lTl lower 
left), motion occurs mostly through ID channels that percolate through the sample. 
Importantly, the channel structures change very slowly with time: some channels close 
and others open up. Since transport is dominated by a small number of channels, this 
give rise to large fluctuations in the conductivity and strong correlations. Therefore, 
in this model, the large noise is due to dynamical inhomogeneities, similar to obser- 
vations in other glasses ( [Glotzer, 2000| [Richert, 2002[ ) . At the lowest , the channels 
disappear and transport occurs via jumps between localized regions. This corresponds 
to the deep insulating regime. In addition, a completely different technique, namely, 
a Monte Carlo simulation of a 3D Coulomb glass out of equilibrium has also found 



evidence for similar dynamical heterogeneities (Kolton et ai, 2005). 

There have been few experimental attempts so far to probe the charge dynamics 
near the MIT in other materials. The most exciting and important has been the study 
of the resistance noise in bulk doped Si ( |Kar et ai, 2003[ ) , which has been a prototypi- 



cal system for investigating the critical behavior near the MIT (Miranda and Dobrosavljevic, 2005 



Sarachik, 1995[ ) for several decades. The results are strikingly similar to those on a 



2DES in Si, namely: i) the noise power increases by several orders of magnitude at ric 
[Fig. lO-lSf a)]. ii) as T decreases, the noise magnitude increases essentially exponen- 
tially for dopings n/uc < 1 (Fig. lO.lSl a) right inset), in) near the MIT, the exponent 
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Fig. 0.18 Resistance noise in P-doped Si. Adapted from Kar et al. (20031. (a) Noise magni- 
tude as a function of doping at different T. Left inset: the same grapii, magnified around ric. 
Right inset: the temperature dependence of the noise magnitude vs. doping, (b) Exponent 
a vs. doping. Inset: the power spectrum S{f) oc 1//". (c) Exponent (1 — /3), a measure of 
correlations, vs. doping. 



a rapidly increases to a value much larger than 1 [Fig. lO.lST b)]. and iv) the exponent 
(1 — /3) of the second spectrum becomes strongly nonwhite (7^ 0) near ric [Fig. IO.18T c)]. 
indicating the onset of correlated dynamics. The similarity of the non-Gaussian spec- 
tra observed in both 2DES and in 3D doped Si near the MIT supports the intriguing 
possibility that such correlated dynamics may indeed be a universal feature of the 
MIT regardless of the dimensionality. 

Resistance noise was also measured in very clean 2D hole systems (2DHS) in GaAs 



(Leturcq et al., 2003 Deville et al., 2005 Deville et al., 20061. Even though an orders 



of magnitude enhancement of the noise was observed with decreasing T and hole den- 
sity Ps, similar to the results on 2DES in Si and bulk doped Si, no evidence for strong 
correlations or glassiness was found. It is likely that the (lack of) disorder here plays a 
crucial role, since glassiness is not expected theoretically in sufficiently clean systems 
( Dobrosavljevic et al., 2003] ) (Fig. I0.16p . It is interesting, though, that the behavior 
of transport and noise in an intermediate regime of ps was attributed to the coex- 
istence of nanoscale conducting and insulating regions, consistent with other studies 
that favor the percolation picture in 2DIIS and 2DES in GaAs piani et al., 2000] 

Allison et al., 200^. The local probes provide ev- 



Ilani et al, 2001 



Gao et al., 2002 



idence (Hani et al, 2000 Hani et al., 2001) that as the density approaches the criti 



cal hole density Pc from the metallic side, the 2DIIS fragments into localized charge 
configurations that are distributed in space, so that the insulating phase is spatially 
inhomogeneous. The complicated structure emerges already on the metallic side of the 
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MIT, reminiscent of the onset of glassiness in disordered 2DES in Si. Clearly, more 
work is needed to examine the charge dynamics in these systems, especially in the 
presence of a higher amount of disorder. 

The resistance noise techniques, long used in studies of various glassy systems, 
are starting to be recognized as powerful tools also in the investigations of the dy- 
namics of cuprates ( Raicevic et al., 2007| [Raicevic et al, 2008[ [Raicevic et ai, 2011[ 
[Bonetti et ai, 2004| Fruchter et al., 2007[ [Caplan et al, 2010[ ). The currently avail- 



able data, however, are scarce and difficult to compare because they probe very differ- 
ent materials and regions of phase space. This warrants further systematic studies of 
the charge dynamics in various cuprates as a function of all the relevant parameters, 
such as doping and disorder, as well as different growth techniques. Comparative stud- 
ies of charge dynamics in a 2DES in Si, as an excellent model system for the MIT and 
out-of-equilibrium behavior, and in more complex materials, such as cuprates, should 
provide important information on these two fundamental problems of condensed mat- 
ter physics. 
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